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We present the X-ray catalog and basic results from our Chandra Large Area 
I Synoptic X-ray Survey (CLASXS) of the Lockman Hole-Northwest field. Our 9 

' ACIS-I fields cover a contiguous solid angle of ~ 0.4 deg^ and reach fluxes of 

5x 10~^^ ergscm"^ s~^ (0.4—2 keV) and 3x 10~^^ ergscm"^ s~^ (2 — 8 keV). Our survey 
^ ■ bridges the gap between ultradeep pencil-beam surveys, such as the Chandra Deep Fields 

. (CDFs), and shallower, large area surveys, allowing a better probe of the X-ray sources 

, that contribute most of the 2-10 keV cosmic X-ray background (CXB). We find a total of 

I 525 X-ray point sources and 4 extended sources. At ^ 10~^^ ergscm"^ s'^^ (2 — 8 keV), 

I our number counts are significantly higher than those of several non-contiguous, large 

' area surveys. Such a large difference is an indication of clustering in the X-ray sources. 

On the other hand, the integrated flux from the CLASXS field, combined with ASCA 
and Chandra ultradeep surveys, is consistent with results from other large area surveys, 
within the variance of the CXB. 



, We see spectral evolution in the hardening of the sources at fluxes below 

^ ' 10~^^ ergs cm~^ , which agrees with previous observations from Chandra and XMM- 

Newton. About 1/3 of the sources in the CLASXS fleld have multiple observations, 
^ ' allowing variability tests. Above 4 x 10~^^ ergs cm~^ s~"^ (0.4 — 8 keV), ~ 60% of 

' the sources are variable. We also investigated the spectral variability of the variable 

sources. While most show spectral softening with increasing flux, or no signiflcant 
spectral change, there are a few sources that show a different trend. 

Four extended sources in CLASXS is consistent with the previously measured LogN- 
LogS of galaxy clusters. Using X-ray spectra and optical colors, we argue that 3 of the 
4 extended sources are galaxy clusters or galaxy groups. We report the discovery of a 
gravitational lensing arc associated with one of these sources. Using red sequence and 
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brightest cluster galaxy methods, we find that the redshifts of the extended sources are 
in the range z ~ 0.5 — 1. The inferred masses within the Einstein radii are consistent 
with the mass profiles of local groups scaled to the same virial radii. 

Subject headings: cosmology:observations — galaxies: active — X-ray:diffuse background 

1. Introduction 

With several ultradeep Chandra and XMM-Newton surveys reaching flux limits as deep as 
/2-10 keV ~ 10~^^ ergs cm~2 s~^ (Alexander et al. 2003; Giacconi et al. 2002; Hasinger 2003), it 
is now clear that point sources account for almost all of the X-ray radiation above 2 keV in the 
universe (Mushotzky et al. 2000; Moretti et al. 2003). These point sources are believed to be 
primarily active galactic nuclei (AGN), based on their X-ray luminosities. It has also been found 
that only ~ 30% of hard X-ray selected AGN show broad lines in the optical, while > 50% appear 
to be normal galaxies in the optical at the sensitivity limits of the current optical spectroscopic 
follow-up (e.g., Barger et al. 2001, 2003). In contrast, the soft X-ray selected AGN from ROSAT 
surveys are mostly identified as broad-line AGN in the optical band (Schmidt et al. 1998). 

The ultradeep surveys cover very small solid angles. In the case of the Chandra Deep Fields 
(CDFs), the combined sky coverage is ~ 0.2 deg^. About 40% variance between fields is seen in 
the integrated fluxes in the 2 — 8 keV band (Cowie et al. 2002), likely as a result of the underlying 
large scale structure. To determine the fractional contribution of point sources to the cosmic 
X-ray background (CXB) with enough accuracy, and to understand how the CXB sources trace 
the large scale structure, a sufficiently large solid angle is needed. While very large area surveys 
exist above 10~^^ ergs cm~^ s~^ (2-8 keV) from ASCA (Akiyama et al. 2003), the data around 
10~^^ ergs cm~^ s~^ , where the point source contribution to the CXB peaks, is limited. 

Several intermediate, wide-field, serendipitous Chandra/XMM-Newton surveys (Baldi et al. 
2002; Kim et al. 2004; Harrison et al. 2003) were designed to increase the solid angle to several 
degrees at a 2 — 8 keV flux limit of 10~^^ ergscm~^s~^ . One of the advantages of such surveys 
is that they sample randomly across the sky, so the probability of all of them hitting overdense 
or underdense regions is small. This is useful in determining the normalization of LogN-LogS. On 
the other hand, serendipitous surveys suffer from the non-uniform observing conditions for each 
pointing, and in most cases, the pointings contain bright sources. The biases introduced by these 
non-uniformities are hard to quantify. The serendipitous surveys also have little power in addressing 
the question of large scale structure traced by X-ray selected AGN, due to the small solid angle of 
each pointing, sparse and random positions on the sky, and the non-uniformity of the observations. 
With serendipitous surveys, it is also difficult to perform extensive optical spectroscopic follow-up 
observations, which are critical in obtaining the redshifts and spectral classifications of the X-ray 
sources. This is due in part to the advent of large format detectors for imaging and spectroscopy 
(like those on the Subaru and Keck telescopes), which are more efficient at targeting large-area. 
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contiguous X-ray surveys, rather than many isolated ACIS-I pointings. 

A contiguous, large solid angle survey can compensate for these disadvantages and bridge 
the gap between the ultradeep "pencil beam" surveys and the large area serendipitous surveys in 
determining both the normalization of the LogN-LogS and the large scale structure. 

In 2001, we began the Chandra Large Area Synoptic X-ray Survey (CLASXS) of the multi- 
wavelength data-rich ISO Lockman Hole-Northwest (LHNW) region. The survey currently covers 
a solid angle of 0.4 deg^ and is sensitive to a factor of 2 — 3 below the "knee" of the 2 — 8 keV 
LogN-LogS. Such a choice of solid angle and depth maximizes the detection efficiency with Chan- 
dra. The large solid angle is important for obtaining statistically significant source counts at the 
"knee" of the LogN-LogS and to test for variance of the number counts on larger solid angles. The 
choice of solid angle is based on the ASCA results that the rms variance of the 2—10 keV CXB 
on a scale of 0.5 deg^ is ~ 6% (Kushino et al. 2002). The expected variance at the angular scale 
of our field should be less than the uncertainty of the CXB flux. The uniform nature of the survey 
allows an unbiased measurement of AGN clustering. We have shown that the ficld-to-ficld variance 
seen between the deep fields can be reproduced within the 9 fields in our survey (Yang et al. 2003; 
hereafter Yang03). 

Our survey region is covered by the deepest 90 and 170/xm ISOPHOT observations (Kawara et 
al. 2004), as well as abundant multiwavelength observations, including the planned Spitzer Space 
Telescope (SST) Wide-Area Infrared Extragalactic Survey (SWIRE). We performed extensive op- 
tical follow-up observations using Subaru, CFHT, WIYN, and Keck to obtain multicolor images 
and spectra of the X-ray sources (Stcffcn ct al. 2004, hereafter Steffen04). These observations 
provide critical information on the redshifts, spectroscopic classifications and luminosities of the 
X-ray sources, as well as on the morphologies of the host galaxies. We emphasize that the subarc- 
second spatial resolution of Chandra is critical to allow the unambiguous identification of optical 
counterparts. 

We describe our observations and data analysis methods in §2 and present our X-ray catalog 
in §3. In §4, we present the X-ray spectral properties and variability of the point sources. In 
§5, we present our analysis of the extended sources and report on the discovery of a gravitational 
lensing arc associated with one of the clusters. We summarize the paper in §6. A companion 
paper by Steffen04 presents our multiwavelength observations and analysis. We will present the 
spatial correlation functions of the X-ray sources in a subsequent paper. Throughout this paper, 
we assume Hq = 71 and a fiat universe with VLm = 0.27 and J^a = 0.73. 
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2. Observations and Data Reduction 

2.1. X-ray Observations 

We surveyed the LHNW field centered at a = 10^34"^, 6 = 57° 40' (J2000). The region has 
the lowest Galactic absorption {Nh = 5.72 x 10^^ cm"^ . Dickey & Lockman 1990). All 9 ACIS-I 
observations were obtained with the standard configuration. The pointings are separated from each 
other by 10' (Figure 1). The fields are labeled LHNWl-9 for reference hereafter. The overlapping of 
the fields allows a uniform sky coverage, because the sensitivity of the telescope drops significantly 
at large ofF-axis angles. Fields LHNWl-3 were observed during April 30th to May 17th 2001, and 
the rest of the fields were observed during April 29th to May 4th 2002. All fields except LHNWl 
have exposure times of ~ 40 ks. LHNWl is located at the center of the field and has an exposure 
time of 73 ks. The observations are summarized in Table 1. 

We reduced the data with CIAO v2.3 and the calibration files in CALDB v2.20. For our 
spectral analysis, we updated the data reduction with CIAO v3.01 and CALDB 2.23 to allow the 
use of CTI corrected calibration files. We followed the CIAO analysis threads (available online 
at http://asc.harvard.edu/ciao/) in reducing the data, including the correction of known aspect 
problems, CTI problems, and removing high background intervals. Background flares were found 
in LHNW3 and LHNW6 and have been removed. The resulting event lists were rebinned into 
0.4 — 2 keV (soft), 2 — 8 keV (hard), and 0.4 — 8 keV (full) broadband images. Spectral weighted 
exposure maps were made for each band for each observation, using the observation specified bad 
pixel files. 

2.2. Source detection 

The detection sensitivity of Chandra drops rapidly beyond 6' off-axis. For this reason, we 
overlapped our ACIS-I fields so that the sensitivity of the survey would be uniform across the 
field. Since the added signal-to-noise from merging the observations is relatively small, we chose 
to detect sources in each observation individually and merge the catalogs, rather than to detect 
sources directly on the merged image. This method certainly loses some sensitivity for very dim 
sources at some locations. However, since our major interest is to obtain a uniform sample for 
statistical and follow-up purposes, such a choice is justified. The method also simplifies the source 
flux extraction because the PSF information could easily be used. Multiple detections of sources 
in independent observations are very useful for checking and improving the X-ray positions of the 
sources. Multiple detections also provide an opportunity for measuring the variability of these 
sources. 

From the various detection tools available, we chose to use the wavdetect tool included in the 
CIAO package (Freeman et al. 2002). Since wavdetect uses a set of scales to optimize the source 
detection, the tool is excellent in separating nearby sources in crowded fields. In general, the 
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method also provides better sensitivity than the classical "sliding box" methods. The drawbacks 
of wavdetect are that it runs rather slowly on large images like the full resolution ACIS-I images, 
and it requires fine-tuning of the parameters. 

We ran wavdetect on the full resolution images with wavelet scales of 1, V^, 2, 2\/2, 4, 4-\/2, 8. 
Although using larger scale sizes could help to detect very far off- axis sources, it is not very useful 
for our survey, because of the overlapping of fields. It also increases the computation time to use 
a large number of scales. We chose to use a significance threshold of 10~^, which translates to 
a probability of false detection of 0.4 per ACIS-I field based on Monte Carlo simulation results 
(Freeman et al. 2002). 

2.3. Source positions 

Observations performed before 2002-May-02 suffer from an systematic aspect offset as large 
as 2" from an error in the pipeline software. This systematic error was carefully calibrated by the 
Chandra X-ray center and corrections are provided. For the affected fields, LHNWl, 2, 3, 4, 5, and 
7, we corrected this error following the standard procedures (see Chandra analysis thread online at 
http : / / asc. harvard .edu/ciao/threads/ arcsec_ correction / ) . 

We further matched the small off-axis X-ray positions reported in each field from wavdetect to 
the optical images. Corrections were then found to maximize the matches. Such corrections are 
very small. The astrometric improvement also only marginally improved the matching between the 
X-ray catalogs, thanks to the excellent astrometric accuracy of the instrument. The corrected X-ray 
catalogs from each observation were then merged (§3.1). A further absolute astrometric correction 
was applied to the merged X-ray catalog to match to the radio sources in the field. In Figure 2, we 
show the offsets between the corrected positions and the optical positions for sources in different 
off-axis angle ranges. Large offsets often occur when the X-ray source is at large oflF-axis angles. 
Offsets are generally small, with a dispersion of 0.23" for all the sources detected. 

2.4. Source fluxes 

Wavdetect is excellent at detecting sources, but it is not always the best method for flux ex- 
traction. Three issues could contribute to an incorrect estimation of source counts in wavdetect. 
First, the flux measurements in wavdetect use a monochromatic PSF size, which, by default, cor- 
responds to an enclosed energy of 0.393 at the energy of choice, or the la integrated volume of a 
normalized two-dimensional Gaussian. Though this parameter is adjustable, larger enclosed energy 
values could cause confusion of close sources. Since the construction of source cells is carried out by 
convolving the source image with wavelet functions, the "smearing" effects of the convolution can 
in general make the source cell large enough to include most of the source photons, but the fraction 
of the flux recovered varies from source to source. Second, due to the statistical fluctuations in 
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the source photon distribution, some sources show multiple peeiks in the convolved image. Unless 
perfect PSF information is available, randomness should exist in determining which peak belongs 
to a single source. This problem is particularly severe when the source is very off-axis and the PSF 
shape cannot be approximated by a Gaussian. The third issue is the background determination, 
a problem other methods also share. The background in wavdetect is obtained in the immediate 
neighborhood of the source. This is useful because of the known large background fluctuations. 
However, if the background is drawn too close to the source, the PSF wing would likely be taken as 
background. This could result in an over-subtraction of the background and lead to underestimated 
source counts. This effect is seen in a correlation of source counts with background density in the 
wavdetect results. All of these issues would lead to an underestimation of source counts. This 
has been noticed in the analysis of the deep Chandra fields (Giacconi et al. 2002; Hornschemeier, 
private communication). 

Because of the spectral diflFerences of the sources and the sensitivity differences between energy 
bands, sources detected in one band are not always detected in another at high significance. There 
is no simple way within wavdetect to provide upper limits for these sources. To obtain the source 
fluxes or upper limits in the non-detection band, an alternative flux extraction method is needed. 

For these reasons, we wrote an aperture photometry tool for source flux extraction. The 
method uses a simple circular aperture which matches the size of the PSF. To do this, we first 
compared the broadband PSFs derived from our observations with the PSF size file provided with 
CIAO, as described below. 

24.1. Broadband PSFs 

Both the PSF library used by the CIAO tool mkpsf and the circularly averged PSFs used by the 
detection codes (psfsize20010416.fits) are generated at monochromatic energies using simulations 
of the telescope. Spectral weighted average energy is usually used for selecting the PSF file for 
broadband images. Since the spectra of the X-ray sources are mostly unknown, an average spectrum 
has to be assumed. Whether such selected PSFs agree with the observed broadband PSFs needs 
to be tested. We constructed "average PSFs" for different off-axis angles using sources which have 
no neighbors within 40" in our 9 observations (Figure 3). It should be noted that these PSFs 
are inaccurate at large scales because the PSF wings, which span more than 1', could not be well 
determined in these observations. The source images from the same off-axis annuli are stacked, 
and the curves-of-growth are constructed. The background regions arc fitted with quadratic forms 
using nonlinear least-square fits. To compare with the library PSFs used by wavdetect, we linearly 
interpolated the library PSFs to the off-axis angles and the spectral weighted averaged energies. 
To account for the fact that part of the PSF wings had been fitted as background in our data, 
we did the same "background fitting" on the interpolated PSFs. This allows a comparison of the 
observed curve with the interpolated PSF. The broadband PSFs are generally narrower than the 
interpolated PSFs, except for one case in the hard band where the off-axis angle is large. 
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2.4-2. Aperture photometry 

We performed our flux extractions in the following way. We used circular extraction cells, 
choosing the radius of cells from the PSF library at a nominal enclosed energy of ~ 95% (the true 
enclosed energy should be >95% based on the discussion above) if the cell size is > 2.5". For 
source close to the aim point, a fixed 2.5" radius was used. The background is estimated in an 
annulus region with an area 4 times as big as the source cell area, with inner radius 5" larger 
than the source cell radius. To avoid nearby sources being included in the background region, 
the background region is divided into 8 equal-sized segments (Figure 4). The mean background 
counts are estimated, excluding the segment which contains the highest number of events. Then 
the 3(7 Poisson upper limit is derived using the approximations provided in Gehrels (1986). The 
background is then recalcTilatcd with only the background segments that contain counts less than 
the upper limit. The net counts are obtained by subtracting the background from the source 
counts within the source cell. We compare the obtained net counts with the net counts obtained 
with wavdetect (Figure 5). While they mostly agree, the source photons derived from our method 
are, on average, higher than those from wavdetect, especially for low-count sources. The average 
increases are 4%, 7%, and 8% for the soft, full, and hard bands, respectively. We hand-checked 
the sources with large discrepancies from the two methods, and we found our estimates to be more 
reliable. 

2.4.3. Exposure time and flux conversion 

The prerequisite for using exposure maps is that the effective area is only weakly dependent on 
energy. This is not the case for our broadband images, where the effective area changes rapidly with 
energy. Using exposure maps blindly, even the spectrally weighted ones, will inevitably introduce 
large errors in the resulting fluxes. However, the vignetting (the positional changes of sensitivity) 
is less sensitive to energy. In other words, if we normalize the exposure maps obtained at different 
energies to the aim points, then the differences between such "normalized exposure maps" are very 
small. 

Based on this fact, we use the exposure maps only to correct for vignetting and compute the 
flux conversion at the aim point using spectral modeling. We first make full resolution spectrally 
weighted exposure maps (using monochromatic maps do not change the results significantly). For 
each source, the exposure map is convolved with the PSF generated using mkpsf and normalized 
to the exposure time at the aim point. This is the effective exposure time if the source is at the 
aim point. 

The conversion factor is then obtained at the aim point by assuming the source has a galactic, 
absorbed, single power-law spectrum. The power-law index is calculated using the hardness ratio 

of each source, defined as HR = Chard/ Cgofti where Chard and Csoft are the count rates in the hard 
and soft bands. XSPEC was used in computing the conversion from HR to F and for calculating 
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the conversions. The degradation of quantum efficiency during the flight of the observatory has 
been accounted for using the standard procedure. 

3. Catalog 

3.1. Merging Catalogs 

We first merged the three band catalogs. We used a 3cr error ellipse from the wavdetect output 
as the identification cell. Flux extraction was then performed on all entries in the merged catalogs 
in all bands using the best position of the sources. We compared the three band catalogs with the 
optical catalog to find the astrometric corrections for each observation, as described in §2.3. The 
nine catalogs were then merged. The fluxes of the sources with more than one detection in the 
9 fields were taken from the observation in which the effective area of the source was the largest, 
except for those sources with more than 2 detections having normalized areas > 80%, where we took 
the averaged flux. We visually checked the final catalog to ensure the correctness of the merging 
process. The final catalog contains 525 sources. 

The distribution of the source off-axis angles in the merged catalog is shown in Figure 6. It 
can be seen that most of the sources fall within the < 6' range. Figure 7 shows the distribution of 
sources with multiple detections. About 1/3 of the sources have more than one observation. 

3.2. Column-by-column description of the catalog 

We present the final catalog in two tables (Table 2a and Table 2b). In Table 2a, we list the 
source positions, fiuxes, and hardness ratios. In Table 2b, we list the source net counts, effective 
exposures, and detection information. 

Table 2a: Basic properties 

Column 1: Source number used in the catalog. The numbers correspond to ascending order 
of right ascension. 

Column 2: Source name follows the lAU convention and should read CXCCLASXS, plus the 
name given in the table. 

Columns 3 — 4: The X-ray position, corrected for the aspect errors of the telescope, if 

applicable, and for the general astrometric solution by comparing with the optical images (see 
§2.3). For sources with multiple detections in the three bands and the 9 observations, the best 
position is taken. 

Columns 5 — 6: Statistical error of the X-ray position quoted from the wavdetect lists. 



-9- 



Columns 7 — 9 : X-ray fluxes in the soft, hard, and fuU bands in units of 10~^^ ergs cm~^ . 

If a source is detected in multiple observations, and if there are more than one observation in which 
the source effective area is more than 80% of the efi^ective area at the aim point, then the mean flux 
is used; otherwise, the flux from the observation that has the largest effective area is used. The 
errors quoted are the la upper and lower limits, using the approximations from Gehrels (1986). For 
any source detected in one band but with a very weak signal in another, the background subtracted 
flux could be negative. In this case, only the upper limit is quoted. 

Columns 10: Hardness ratio as defined in §2.4.3. The upper or lower limit is listed for a 
source with no net extracted photons in the soft or hard bands. 

Table 2b: Additional Properties 

Column 1: Source number. 

Columns 2—4: Net counts in the soft, hard, and full bands. If a source is detected in multiple 
observations, then the observation in which the source has the largest effective area (see Column 
8) is used. As in Table 2a, for the sources with negative counts, only the upper limits are listed. 

Columns 5—7: Effective exposure time in each of the three energy bands from the exposure 
map. 

Columns 8—9: Detection information. Column 8 is the LHNW field number where the source 
has the largest effective area. Column 9 lists the LHNW field numbers (each digit represents a field 
number) in which the source has been detected in at least one of the three bands. Sources with 
multiple detections are necessary for the detection of variability (see §5.2). 

4. Results 
4.1. Number counts 

4.1.1. Incompleteness and Eddington Bias 

Incompleteness can be caused by energy or positional dependence of the sensitivity of X-ray 
telescopes. Because the spectrum of a source carries important information on the physical nature 
of the source itself, sources of different spectra are usually categorized as different types. The energy 
dependent sensitivity acts like a filter in selecting "hard" and "soft" types of X-ray sources. The 
soft band detected sources always contain more soft spectrum objects than the hard band detected 
sources and vice versa. Unless the fraction of each type remains constant for all fluxes (which we 
now know is not true), the energy dependent incompleteness cannot be easily corrected. This issue 
is very important in interpreting the fraction of different types of objects in flux limited surveys. It 
is desirable to obtain number counts for each type of source, but it is hard to do that for the CXB 
sources, because the spectra are hard to determine. For our medium deep survey, it is sensible to 
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foUow tradition and only discuss the number counts in energy bands. 

The positional dependent incompleteness is caused by vignetting and aberration of the X-ray 
optics. The vignetting causes the effective area to drop with off-axis angle, and the aberration 
makes the off-axis PSF larger so that it includes a larger number of background events in the 
source cell. The net effect is that the sensitivity of source detection drops with increasing off-axis 
angle. The sky area is therefore flux dependent. 

These effects can be investigated via Monte Carlo simulations. Wc first generated background 
images using observations of fields ^1 and #4, which represent the 70 ks and 40 ks exposures. 
Point sources are removed from the images, and the holes left in the images are filled by sampling 
the local background. Random sources are generated uniformly on the background images. The 
fluxes of the sources are generated by randomly sampling a complete subset of the combined 
Chandra Deep Fields catalog (Alexander et al. 2003). The subset contains only sources with hard 
band fluxes > 5 x 10~^^ ergs cm~^ s~^ and effective exposure times > 200 ks. The input fluxes 
are converted to on-axis counts assuming power-law spectra with F = 1.4. The exposure map for 
each image is consulted to find the vignetting effect at the source location, and the normalized 
exposure is multiplied by the true counts to obtain the "observed" net counts. Only sources with 
more than 3 counts are used in the simulation to avoid adding too many undetectable dim sources 
to the background. We use the CIAO tool mkpsf to generate realistic source shapes at the source 
positions and energies. The PSF is then sampled to have the same number of photons as in the 
source. We chose to use mkpsf instead of using the Chandra simulator MARX because we find the 
PSF library used by mkpsf better resembles sources at large off-axis angles. The number density 
of the sources is chosen to be 2 times higher than the observed density to increase the number of 
simulated sources without affecting detections. We ran 100 simulations on the 2 fields and the three 
energy bands and detected the sources using wavdetect with identical parameter settings to those 
we used in preparing the observed catalog. Because of the large computation time, the number of 
simulations that can be done is limited. 

We then compared the output catalogs with the input source catalogs. Because of the small 
size of the simulation, the completeness within 4" is not well determined, and a 5% uncertainty 
exists in the determined fractions. Fortunately, the PSF effect is small at such small off-axis angles. 
For a given flux threshold, the fraction of source detections drops monotonically with off-axis angle. 
This relation is fitted between 4' and 10' with a linear least-squares fit. The 95% complete off-axis 
angle limit is then taken from the interpolation of the fit. The resulting flux thresholds versus 
off-axis radii arc shown for the three bands and the two exposure times in Figure 8. We note that 
at large off-axis angles, the sensitivity drops rapidly. This is due to the choice of wavelet scales. 
When the largest scale used becomes smaller than the PSF size of the source, wavdetect is no longer 
sensitive. This effect, however, is not important for our observations, because most of the sources 
of interest are within 6' off-axis, thanks to the overlapping of fields. With these curves, we are 
able to make threshold maps for the combined catalogs of all three bands. Sources at very large 
off-axis angles are excluded from the study of the LogN-LogS. The combined solid angle versus flux 
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thresholds is shown in Figure 9. 

The Poisson fluctuations in the source fluxes could result in an overestimation of number counts 
close to the detection limits. This is known as the Eddington Bias. The effect depends both on the 
slope of the LogN-LogS and the level of fluctuation. For the CLASXS fleld, the detection threshold 
is below the "knee" of the LogN-LogS, and the Eddington bias is relatively small. We corrected this 
bias using the method described in Viklinin et al. (1995). In Figure 10, we compare the average 
input flux with the average output flux at different off-axis angles from the simulations. For the 
soft band, the correction is only important below 2 x 10~^^ ergscm~^s~^ . For the hard band, the 
correction is important below 8 x 10~^^ ergscm~^s~^ . We fit flux-flux curves in Figure 10 for the 
different off-axis angles with fourth order polynomials and correct the source fluxes in the observed 
catalog using these fits. 



4.1.2. Number counts 

Sources are selected by consulting the threshold map at the source positions and including 
only those with Eddington bias corrected fluxes higher than the threshold map values. Sources 
very far off-axis are excluded from the analysis. With these selections, we used a total of 310 and 
235 sources in the soft and hard bands, respectively, to construct the LogN-LogS. The cumulative 
LogN-LogS relations are computed using the formula 

where O is the complete solid angle. We show the results in Figure 11 in the soft and hard bands 
with 1(7 Poisson errors. The differential LogN-LogS for the two bands are shown in Figure 12, 
which are calculated using the formula 

— =y^^ (2) 

in units of deg^^ per 10~^^ ergscm~^s~^. We fit the resulting differential number counts with 
single or broken power laws in the form of 



dN 

Is ""'10 



^ = ^o(t7^)-" (3) 



using error weighted least-square fits. Since our survey best samples the "knee" of the LogN-LogS, 
the slope of the power-laws are not well constrained due to the lack of data points both far above 
and below the "knee". On the other hand, Uq is better determined, to within 1%. 

For the soft band, we fit the number counts between 10~^^ and 10~^^ ergs cm~^ s~^ with a 
power law. We find the best-fit parameters to be o; = 1.7 it 0.2 and no = 12.49 it 0.02. The slope 
is in good agreement with previous observations, such as the Chandra Deep Field-North (1.6 it 0.1, 
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Brandt et al. 2001), SSA13 (1.7 ± 0.2, Mushotzky et al. 2000), and the compiled wide fields from 
Chandra, XMM-Newton, ROSAT, and ASCA {imto'^l, Moretti et al. 2003; hereafter, Morctti03). 
The normalization also shows excellent agreement with the compiled results from the large area 
survey of MorettiOS, which has an effective solid angle at 10~^^ ergs cm~^ s~^ larger than that 
of CLASXS. Above 10~^^ ergscm~^s~^ , the slope steepens, but the fluctuations in the number 
counts make it difficult to find a reasonable fit. However, the LogN-LogS is apparently consistent 
with a slope of a = 2.5, shown as the dotted line at these fiuxes. 

Similarly, we model the hard band number counts with a broken power law and obtain the 
following best-fit parameters. For S > 10~^^ ergscm"^ , a = 2.4 it 0.6 and no = 45.6 ± 0.5; for 
3x 10~^^ < S < 2x 10~^^ ergscm"^ s~^ , a = 1.65±0.4 and no = 38.1±0.2. For comparison, we also 
plot the best-fit cumulative LogN-LogS from Moretti03 and differential LogN-LogS from the SEXSI 
fields (Harrison et al. 2003) and from Cowie et al. (2002). At fluxes below 8 x 10~^^ ergscm"^ s~^ , 
the differential LogN-LogS for all the flelds agrees within the errors. The difference in the total 
counts at a flux limit between the CLASXS fleld and the MorettiOS flelds is also small. An apparent 
difference is seen around 10~^^ ergscm~^s~^ : the total counts at 10~^^ ergs cm~^ s~^ are ~ 70% 
higher than those from Moretti03. This is significant at greater than the 3cr level. 

The integrated flux between 3 x 10~^^ and 8 x 10~^^ ergs cm~^ s~^ from the LogN-LogS is 
(1.2 ± 0.1) X 10~^^ ergs cm~^ s~^ deg~^. This is ~ 20% higher than that from Moretti03 and 
SEXSI in the same flux range. Since there is little difference in the number counts between the 
CLASXS fields and the other large solid angle surveys at fluxes lower than 8 x 10~^^ ergs cm^^ s~^ , 
we should expect little difference below the survey limit on the same angular scales. If integrated to 
lower fluxes, and including the integration from ASCA above 8 x 10~^^ ergs cm~^ , the fractional 
difference between the CLASXS field and the other large solid angle surveys can be reduced to 
~ 10% without considering the possible biases. This difference is higher than expected from the 
variance in the CXB from ASCA observations but is consistent with recent observations with 
RXTE/PCA, where a 7% variance is seen among several ~ 1 deg^ fields (Revnitsev et al. 2003). 

The uncertainty of the hard CXB itself is 10 — 15%. The differences in the integrated point 
source fluxes from various large flelds are within this uncertainty. In terms of the true contribution 
from point sources to the CXB, a field with a solid angle of ~ 0.3 deg^ seems to be large enough 
to be representative. 

The detection of the large variance in the hard band source counts around 10^^*^ ergs cm^^ s^^ on 
our survey scale seems to indicate that the sources that emerge at this flux are more clustered on 
the sky than the soft band selected sources. These sources could account for most of the cosmic 
variance observed. This issue will be discussed in a separate paper. 
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4.2. Spectral properties and variability of the CXB sources 

4.2.1. Hardness ratio 

We employ a hardness ratio (§2.4.3) to quantify statistically the spectra of the CXB sources 
in our field. Figure 13 shows the distribution of hardness ratio versus full band flux. We have also 
marked the hypothetical photon indices (F), assuming the hardness ratio change is purely due to 
the slope change of a single power-law spectrum. At fluxes > 3 x 10~^^ ergs cm~^ s~^ , most sources 
cluster around F ~ 1.7. At lower fluxes, the hardness ratio distribution scatter increases and the 
relative number of hard sources increases. Below 3 x 10~^^ ergs cm~^ s"^ , the data show a paucity 
of hard sources. This is a selection effect caused by the sensitivity in the hard band being lower than 
in the soft band for most spectra. We stacked the sources in flux bins and calculated the hardness 
ratios of the stacked spectra. Figure 14 shows the stacked hardness ratios from both the CLASXS 
> 10~^^ ergs cm~^ s^^ sample and the combined CDFs > 10"^^ ergs cm~^ s~^ sample. The flux 
thresholds are chosen to avoid selection effects caused by the sensitivity differences between the 
soft and hard bands. It is apparent that the results from our data and those from the CDFs agree 
well. 

The spectral flattening at low fluxes has been observed by several authors (e.g., Mushotzky 

et al. 2000; Tozzi ct al. 2001; Piconcelli ct al. 2003; Alexander et al. 2003) with observations of 
different depths. Spectral analyses with XMM-Newton observations indicate that such a flattening 
is mainly caused by absorption. These obscured AGN must dominate the population around the 
"knee" of the LogN-LogS to account for the flat spectrum of the CXB. Since most of the XMM- 
Newton spectral observations have reached a few times 10"^^ ergscm"^ (Piconcelli et al. 2002), 
and the mean spectrum at this threshold is still too soft compared with that of the CXB, a sharp 
increase of obscuration or a change of spectral shape at a flux ^ 10~^^ ergs cm~^ is inevitable. 
Such a sharp change is seen in the change of hardness ratio in our wide-field sample. 

4.2.2. Variability 

X-ray variability is an important factor in distinguishing AGN from starburst galaxies. Almost 
all AGN vary in X-rays, except those sources which are Compton thick. Alexander et al. (2001) 
showed that only a small fraction of the optically faint X-ray sources vary. Possible explanations 
could be that a large fraction of the optically faint sources are Compton thick, or that the amplitude 
of variation of the optically faint sources is much lower than that of the broad and/or narrow- line 
AGN at the same flux thresholds. 

We examine the variability of sources that have been detected in more than one of our obser- 
vations. Since the observations were taken in two groups, separated by about one year, and each 
group of observations were taken within a few days (see Table 1), we are able to test variability on 
timescales of days and/or one year, depending on the location of the source. 



For timing analyses with low counts per bin, the usual statistic is inadequate. We use 
the C-statistic (Cash 1979) in testing the significance of variability. Cash (1979) showed that the 
C-statistic (a reduced form of likelihood ratio) written as 



N 



AC = -2^[niln{ei) 



-Ci- niln{ni) + rii] 



(4) 



is asymptotic to a distribution with A'' — 1 degrees of freedom, where is the observed counts 
in the ith sample, Cj is the expected counts in that sample, and N is the total number of samples 
used. We restricted the sample for the variability test to sources with expected counts greater than 
10 in all observations. The null hypothesis rejection probability was chosen to be 0.01. 

A total of 168 sources were tested for variability, of which 42 sources are significantly variable 
and 28 sources show variability on timescales of days. There are 29, 16, and 30 variable sources de- 
tected in the soft, hard, and full bands, respectively. Figure 15 shows the light-curves of the sources 
that were tested to be variable in any of the three energy bands. In the top panel of Figure 16, we 
show the fraction of variable sources detected versus flux. Between 4 — 8 x 10~^^ ergscm~^s~^ , 
70% of the sources tested show variability. This fraction drops dramatically as the flux decreases 
and, at 10~^^ ergscm^'^s"^ , reaches below 20%. This is at least in part due to the selection ef- 
fect that larger variability is needed at lower fluxes to make the test significant. At fluxes above 
8 X 10~^^ ergscm"^ s~^ , only one of the four sources tested (25%) was found to be variable. 

Following Nandra et al. (1997), we define the magnitude of variability as the "excess variance", 
the error subtracted rms variance 



where fi is the flux in each observation, (i is the mean of the fluxes, and cTj is the Poisson error of 
the flux. By assuming the same power density spectrum of X-ray variability for all AGN, cr^^^ can 
be used as a good indicator of whether the variability exceeds the Poisson noise. It has been found 
that there exists a good anti-correlation between cr^^g and AGN luminosity (Nandra et al. 1997) 
in local AGN samples. 

In Figure 17, we show the excess variance of sources that had been detected to be variable 

versus X-ray flux in the three energy bands. At high fluxes, the average o'^„js is significantly lower 
than at lower fluxes. As mentioned above, variability is harder to detect for low flux sources, unless 
the source is more variable than that of the brighter sources, so this bias could explain why there 
are very few low flux, low variability sources in the plot. In addition, the sources we detect to be 
variable are generally soft. This is consistent with the observation from the CDFs that optically 
faint sources (most of which are hard spectrum AGN) are less variable (Alexander et al. 2001). 




(5) 
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4-2.3. Spectral variability 

Very little is known about the spectral variability of the sources that contribute the most 
to the CXB due to a lack of data. Spectral variability is seen in about half of the well-studied 
brighter sources, with a general trend of softening of the 2—10 keV spectra with increasing source 
intensity. But a counterexample is NGC7469, where the spectrum flattens when the source flux 
increases (Barr 1986). The variability could be accounted for either with a change in the relative 
normalization of the different spectral components or by variation in the absorption. 

In FigTirc 18, we show hardness ratios versus full band fluxes for the variable sources. While 
most of the sources show either no clear spectral variability, or a trend of spectral softening with 
increasing flux, there are a number of sources that clearly become harder with increasing flux. 
There are also a few sources that exhibit a mixed trend. On average, these sources tend to have 
softer spectra with increasing flux. 

5. Extended Sources 
5.1. Detection 

We searched the 0.4—2 keV images of each observation for extended sources using the vtpdetect 
tool provided in CIAO. The method uses Voronoi tessellation and percolation to identify dense 
regions above Poisson noise. This method performs best on smooth overdense regions but could 
confuse crowded point sources. We chose to use a threshold scale factor of 0.8 and a maximum 
probability of false detection of 10~^ and to restrict the number of events per source to > 30. 
We used default values for the rest of the parameters. This choice of parameters maximizes the 
detection of low surface brightness sources at high significance. We visually examined the source 
list to screen out apparent blended point sources. The candidates were then selected by comparing 
the 99% PSF radius with the equivalent radius of the source region, and only sources with a PSF 
ratio (defined as \/\A/'k) /r^^., where A is the area of the source region reported by vtpdetect., and 
rgg is the 99% PSF radius at the off-axis angle) higher than 10 were considered extended (Table 3). 

Four sources were found to be significantly extended, and all but Source 3 have an off-axis 
angle of < 5' in the X-ray observations. Source 3 is at an off-axis angle of 8.4'. With the X- 
ray image alone, one could not rule out the source being a blend of point sources. However, a 
bright gravitational lensing arc found in the optical image (see §5.6) at the X-ray peak makes it 
very likely that the X-ray emission is associated with a cluster. Considering the non- uniformity 
of the detection due to vignetting and PSF effects, the number counts for extended sources above 
3.7 X 10~^^ ergs cm~^ s""-*^ are roughly > 10 deg~^. This agrees with the LogN-LogS of clusters at 
these fluxes found in the CDFs (Bauer et al. 2002). It is interesting to note that all 4 extended 
sources are found on only two of the overlapping ACIS-I fields in the north of the LHNW region. 
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5.2. Optical imaging observations 

Optical observations arc describe in detail in Steffen04. Optical images were obtained with 
Suprime-Cam on Subaru (in B, V, R, I, z' ) and with the CFH12K camera on the Canada-Prance- 
Hawaii Telescope (CFHT; in R, B, and CFHT z'). The 2a limits in B, V, R, I, and z' are 27.8, 
27.5, 27.9, 26.4, and 26.2. 

X-ray contours overlaid on R band optical images are shown in Figure 19. We examined the 
number counts of galaxies within circular cells with fixed radii of 0.5'. At a threshold of i? < 24, a 

total of 19, 28, 18, and 9 galaxies were found within the cells centered at the X-ray peaks of each 
extended source. Because a star is found at 0.6' south-east of the X-ray peak of Source 2, the 
galaxy counts could be underestimated. Compared with the expected 6.7 galaxies per cell obtained 
from the whole field, the over densities of galaxies in Sources 1 and 2 are > 3a, while the over density 
of galaxies in Source 3 is ~ 3cr. Source 4 does not show significant clustering of galaxies in the R 
band image. Sources 1 and 2 are very close to each other, with a separation of ~ 2'. The closeness 
and the elongated morphology of the two sources suggest that they are undergoing interactions. 
Source 3 is extended along the east-west direction with multiple peaks. All 4 sources show bright 
elliptical galaxies at the X-ray peaks. 

5.3. X-ray spectra 

We extract very coarse spectra (grouped to > 15 counts per bin to allow the use of the 
statistic) and attempt to constrain the properties of the clusters. We fit the data with a simple 
MEKAL model in XSPEC (vll.2), with a fixed abundance of 0.3 of the solar value and a fixed 
Galactic absorption. We restrict the spectral fitting to within 0.5 — 5 keV, because the signal-to- 
noise ratio is poor outside of this range. 

The source extraction and background regions of Sources 1 and 2 are shown in Figure 20. 

The regions avoid the point sources between the two clusters. The spectra are shown in Figure 
21. For Source 1, we found the best fit to be kT= 1.4]to2 keV and a = 0. 51^02, with a reduced 
= 8.8 for 9 degrees of freedom. This agrees with the redshift estimates using the optical data 
(Table 4). Fitting the same model to the spectrum of Source 2 with the redshift fixed to z = 0.5 
yields kT= 3.lt}^6^ keV, with a reduced = 4.6 for 9 degrees of freedom. The constraint on 
the temperature is poor, but the probability that the temperature of Source 2 is significantly 
different than that of Source 1 is low. This can be seen in Figure 22, where the joint probability 
contour of the temperature from the two sources is shown. The confidence level for the two sources 
having different temperatures is only 2a. Combining the two data sets and fixing z = 0.5, we find 
kT= l.ltli keV. 

The spectrum of Source 3 shown in Figure 23 was extracted from a circular region with radius 
36". The background was extracted from an annulus with inner radius 36" and outer radius 60". 
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The data cannot constrain the model very well, but a simple fit with an absorbed power-law shows 
that the spectrum is very soft with photon index F = 2.6 and reduced = 7.2 for 7 degrees of 
freedom. Fitting with a MEKAL model and assuming a redshift of z = 1.1 (see §5.5), we obtain 
a temperature of 2.3^'^Q'g with reduced = 0.77. The temperature is insensitive to the redshift 
between z = OA — 1.4. The fact that the MEKAL model fits the data better makes it less likely 
that Source 3 is a blend of several point sources. 

With only 30.7 net counts, it is impossible to model the spectrum for Source 4. However, the 
source has very few counts above 2 keV, indicating that the temperature should be low if the source 
is at z > 0.4, as implied from the optical data. 

The virial masses of the extended sources can be roughly estimated using the best-fit M — L 
relation (Finoguenov et al. 2001), M500 = 2.45 x 10-^^ T^-^^, where M500 is the mass within a radius 
where the overdensity is 500. The results are shown in Table 5. All of the sources belong to low 
mass clusters or groups, and this result is not very sensitive to the redshift because of the very soft 
spectra. 

5.4. Angular sizes 

The angular sizes of the sources were quantified by the widths of the radial profiles. We fitted 
the radial profiles of the sources with integrated 2-D Gaussian curves, which describe the low S /N 
ratio data reasonably well. We constructed the cumulative counts as a function of off-source radius 
(curve-of-growth) . Exposure maps were applied to correct for vignetting. Nearby point sources 
were removed and replaced with background noise. The background regions were selected visually 
and fitted with a quadratic form plus a constant. The curves-of-growth were then normalized to 
the best-fit backgrounds. The normalized curve-of-growth for each source is shown in Figure 24. 
This left only one parameter to be determined — the widths of the curves. The best-fit la radii are 
listed in Table 5. 



5.5. Redshifts 

We infer the redshifts of the extended sources using the red sequence method, as well as 
the brightest cluster galaxy (BCG) method. Based on observations of clusters, there is usually 
a population of early- type galaxies which follow a color-magnitude relation (red sequence). This 
relation changes with redshift in a predictable way, such that a robust two-color photometric redshift 
can be obtained (Gladders & Yee 2000). Color-magnitude plots of the sources within 0.5' of the 
X-ray centers are shown for each extended source in Figure 25. Red sequences can be clearly 
seen in Sources 1, 2, and 4. By comparing with the models from Yee & Gladders (2001), we can 
estimate the redshifts for these three extended sources (Table 4). Source 3 does not show a clear 
red sequence. 
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BCGs are often used as distance indicators, because they have almost constant luminosity 
(Humasom, Mayall, & Sandage 1956). One of the difficulties in applying this method is that with 
optical images alone, it is hard to distinguish between the background and the cluster members, 
unless a density peak can be clearly determined. In our case, this is less worrisome because bright 
spheroidal/lenticular galaxies are found at the X-ray peaks of all of the extended sources. This 
clearly associates these galaxies with the clusters. Furthermore, these galaxies are also the brightest 
early- type galaxies in the regions where X-ray emission is significant. Following Postman &; Lauer 
(1995), we fit the radial profile of each of the BCGs to obtain the magnitudes within angular radius 
rj„ and the slope of the profile {a = dlogL„Jdlogr\j.^, where Lm is the luminosity within Vm. We 
eye examine the profile so that nearby galaxies are not included in the aperture. K correction is 
performed using Kr = 2. 5Zo5io(l + 0.962;). We then solve the redshifts for each galaxy by assuming 
the cosmological parameters in §1. The resulting redshifts are listed in Table 4. 

While it appears that the BCG method produces higher redshifts than the red sequence 
method, the differences are not significant, given the large uncertainties in both methods. The 
redshifts of Sources 1 and 2 also agree with the spectral fitting results from the X-ray data. 

The X-ray luminosities of the extended sources are listed in Table 5, assuming the red sequence- 
determined redshifts except for Source 3, where BCG redshift is adopted. Within errors, the 
temperatures and luminosities of the sources agree with the scaling law found in high-redshift 
X-ray clusters (Ettori et al. 2003), but the constraint is weak. 



5.6. Discovery of a gravitational lensing arc 

We have found a gravitational lensing arc close to Source 3 (Figure 26) . The arc has an angular 
radius of ~ 6" . A bright spheroidal galaxy is clearly associated with the arc. A possible counter 
arc is seen connecting to the west of the bright galaxy but is not fully resolved. With B, V, R, /, 
and z' observations, we can estimate photometric redshifts for the cD galaxy and the arc using the 
publicly available photometric redshift code Hyperz (Bolzonella et al. 2000). We find photometric 
redshifts for the cD galaxy and the arc oi z = 0.45 and z = 1.7, respectively. The redshift of 
the galaxy is slightly different than the redshift of the cluster obtained using the BCG method. 
From our experience, one often needs at least 7 colors to obtain a secure photometric redshift. The 
redshift estimates therefore need verification. We now discuss the estimated lensing properties. 

If the source is at Zgrc ~ 0.45 and the arc is at Zarc ~ 1-7, then we can estimate the mass 
within the Einstein radius (reasonably approximated by the radius of the arc) as 

M{e < Oe) = 1.1 X ioi^(^)2(:^M^)Me , (6) 

where D^, Ds, and D^s are, respectively, the angular diameter distances (in units of Gpc) of 
the lens, source, and the distance between the lens and source. With Zgrc ~ 0.45 and Zarc ~ 1-7, 
we obtain M{9 < 9e) ~ 3.3 x 10^^ M©. We compare this mass with what would be expected if 
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the source were a group of galaxies at z = 0.45, assuming the mass profiles are self-similar. By 
fixing the redshift, the X-ray spectra yield a temperature of 2.2 keV. The virial radius is roughly 
I'bOO = 0.63 X y^kT) = 955 kpc (Finoguenov et al. 2001), where rsoo is defined as the radius 
within which the overdensity is 500. The size of the arc at z ~ 0.45 is TdfQ 37 kpc — 0.036r5oo* 
Comparing with the mass profiles of NGC2563, NGC4325, and NGC2300 (Mushotzky et al. 2003), 
the mass inside the Einstein radius agrees very well with that of a group of galaxies. The virial 
mass of the group can then be estimated to be ~ 1.2 x 10^^ M0 (Finoguenov et al. 2001). 

If the cluster is at 2; ~ 0.7, as implied from the BCG method, and if the best-fit temperature 
kT = 0.23 keV is assumed, then we can search for the best redshift of the lensed galaxy, so that the 
mass within the Einstein radius agrees with the mass profile of groups. We find that if the lensed 
galaxy is at z = 1.8, then the mass within the Einstein radius is M{9 < 6e) ~ 3 x 10^^ Mq, which 
fits the mass profile of groups. 

If the redshift estimate is correct, then the arc system is very similar to the one discovered in 
the ROSAT survey of the Lockman Hole (Hasinger et al. 1998b). High-redshift gravitational 
lensing arcs are rare objects so far observed. However, since our large area survey is very similar 
in sky area and depth to the ROSAT Deep Survey, and since both have produced a detection of a 
strong arc, the probability of detection seems high. Larger area surveys of X-ray selected clusters 
of galaxies with deep optical follow-up would help to determine the probability of detection. Such 
observations should put useful constraints on fi^ and on the density of galaxies at high redshifts 
(Cooray 1999). 

It is interesting to note that all four of our clusters may have redshifts z ^ 0.4 — 0.5 and arc 
located within a region of only ~ 20' at the north-east corner of our field. This corresponds to a 
comoving radius of ~ 5 Mpc. The implications of such large scale structure on the CXB need to 
be investigated further. 

6. Summary 

In this paper, we presented our CLASXS X-ray catalog. Our survey covers a ~ 0.4 deg^ 
contiguous area in an uniform manner and reaches fluxes of 5 x 10~^^ ergs cm""^ in the 0.4—2 keV 
band and 3 x lO"-*^^ ergscm^^ s~^ in the 2 — 8 keV band. We found a total of 525 point sources 
and 4 extended sources. We summarize our results as follows. 

1) The number counts in the 0.4 — 2 keV band agree very well with other large area surveys. On 
the other hand, the number counts in the 2 — 8 keV band deviate significantly from other large area 
surveys at the "knee" of the LogN-LogS, possibly as a result of the underlying large scale structure. 
The total 2 — 8 keV band fiux agrees with the observed CXB fiux within the observed variance of 
the CXB, indicating that the true normalization of the CXB can be determined using fields with 
solid angles ~ 0.3 — 0.4 deg^. 
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2) The hardness ratios of the sources in the CLASXS field show a significant change at /2_8 keV ~ 
10~^^ crgscm^^ s^^ , which bridges the range sampled by previous studies and confirms the results 
found in deep Chandra/XMM-Newton surveys. About 60% of the sources with full band fluxes 
> 4 X 10"^"^ show significant variability, while the fraction drops dramatically with decreasing 
flux, at least partly due to selection eflFects. Most sources show no change of hardness ratio or 
anti-correlation with flux. But some sources show a positive correlation or mixed trends. 

3) We report on the X-ray and multicolor analysis of four extended sources. We argue that the 
sources are likely low mass clusters or groups at redshifts ~ 0.5. Two of the clusters are probably 
interacting or merging. 

4) We report on the discovery of a gravitational lensing arc. The lensing cluster is consistent with 
being at a redshift oi z = 0.45 — 0.7. 
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Tablc 1. Observation Summary 



Target Name 


0:2000 


^2000 


Obs ID 


Sequence # 


Observation date 


Exposure' 


LHNWl 


10 34 00.24 


+57 46 


10.6 


1698 


900057 


05/17/01 


18:29:38 


72.97 ks 


LHNW2 


10 33 19.82 


+57 37 


13.8 


1699 


900058 


04/30/01 


10:59:38 


40.74 ks 


LHNW3 


10 34 36.12 


+57 37 


10.9 


1697 


900056 


05/16/01 


12:46:50 


43.72 ks 


LHNW4 


10 32 04.20 


+57 37 


15.6 


3345 


900184 


04/29/02 


03:23:45 


38.47 ks 


LHNW5 


10 34 00.31 


+57 28 


15.6 


3346 


900185 


04/30/02 


02:03:59 


38.21 ks 


LHNW6 


10 33 20.28 


+57 55 


15.2 


3343 


900182 


05/03/02 


09:11:41 


34.04 ks 


LHNW7 


10 32 44.23 


+57 46 


15.2 


3344 


900183 


05/01/02 


20:03:06 


38.54 ks 


LHNW8 


10 34 36.26 


+57 55 


15.6 


3347 


900186 


05/02/02 


14:16:27 


38.46 ks 


LHNW9 


10 35 14.28 


+57 46 


15.2 


3348 


900187 


05/04/02 


11:01:47 


39.52 ks 



Total good time with dead time correction. 
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Table 3. Extended Sources 



Source # 


02000 




Aa (") 


AS (") 




PSF ratio'' 


Net Counts'^ 


Field"* 


1 


10 35 25.4 


+57 50 48 


2.628 


1.044 


4.786 


37.89 


100.1 ± 11 


9 


2 


10 35 13.4 


+57 50 17 


3.312 


1.476 


4.029 


39.33 


87.7 + 11 


9 


3 


10 35 37.9 


+57 57 15 


3.060 


1.476 


8.422 


19.09 


77.0 ± 10 


8 


4 


10 34 30.8 


+57 59 12 


3.132 


2.196 


4.016 


16.69 


30.7 + 6.6 


8 



^OfF-axis angle in the field the source is detected. 

''Defined as y/{A/'K)/rgg, where A is the area of the source region from the vtpdetect report and rgg 
is the 99% PSF radius at the off-axis angle. 

"^Net counts reported by vtpdetect 

''LHNW field number where the source has the smallest off-axis angle. 

Table 4. Redshift estimates for the extended Sources 



Source # 


zrs 


zbcg 


^X—ray 


1 
2 
3 
4 


0.50 
0.50 

0.45 


U-0O-().08 

0.73ini 
0.451°;°^ 


0.5 + 0.2 

0.5 + 0.2 



Table 5. Properties of the extended sources 



Source it 








core radius (") 


.fo.r, cn;, , \ - 


' L,,„i ■ 


1 
2 
3 
4 


.50 
.50 
.73 
.45 


1 4+0-8 
■■■•^-0.4 
1+6.5 

9 q+1-0 

^•"^-o.g 
1.0 (fixed) 


0.45+°f 

2.0^1^4 

.24 


12.9 
17.0 
14.7 
11.8 


1.6 
1.2 
1.5 
.42 


2.2 
1.5 
5.1 
.45 



''Listed are single parameter la errors, 
''unit: 10" M0 

•^Unit: 10""ergscm-2s"' . 
''Unit: 10*^ ergs s"' . 
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Fig. 1. — Layout of the 9 ACIS-I pointings. Gray scale map shows the adaptively smoothed full 

band image. The exposure maps are added (light gray) to outline the ACIS-I fields. Fields are 
separated by 10' from each other. The field numbers (LHNWl-9) are shown at the center of each 
ACIS-I field. 

Fig. 2. — Offset between the corrected X-ray positions and the optical positions. The 135° and 
45° , as measured from positive x-axis, shaded histogram represents sources with best positions 
within 4' and between 4' and 6' of the optical axis, respectively. The unshaded histogram shows 
the sources with off-axis angles > 6'. 

Fig. 3. — Broadband PSFs obtained from our observations (solid lines) compared with the 
monochromatic PSFs from the PSF library (dashed lines). Both the observed broadband PSFs 
and the monochromatic PSFs are normalized to the wing. From the narrowest to the broadest, 
each broadband PSF is constructed within each of the off-axis angle intervals 0'-4', 4'-6', 6'-7', 
7'-8', 8' 9', 9' 12'. The library PSFs arc taken at the midpoints of these off-axis intervals, (a) Soft 
band PSF vs. 0.91 keV library PSF; (b) hard band PSF vs. 4.2 keV library PSF. 

Fig. 4. — Examples of the source and background regions used in the flux extraction. The smaller 
circle is the source region. The background regions arc shown as segments of an annulus. Segments 
with counts below 3a of the mean are used in the final background estimation and are marked with 
'X' symbols, (a) An isolated source; (b) a source with a close neighbor. 

Fig. 5. — Comparison of net counts from wavdetect and our aperture photometry (marked as 
XPHOTO) for the (a) soft, (b) hard, and (c) full bands. 

Fig. 6. — Distribution of off-axis angles of the best positions. 

Fig. 7. — Distribution of multiple detections. 

Fig. 8. — Thresholds for 95% completeness vs. off-axis angle. Sources with flux and off-axis 
angle combinations above these curves are complete. Solid, dashed, and dotted lines represent 
the threshold curves for the soft, hard, and full bands, respectively. Squares represent the 70 ks 
exposure and diamonds the 40 ks exposure. 

Fig. 9. — Survey effective solid angle vs. flux. Soft band (solid line); full band (dashed line); hard 
band (dotted line). 

Fig. 10. — Average output fluxes from wavdetect vs. average input fluxes for the simulated sources 
at a set of off-axis angles (diamonds). The Eddington bias is seen in the overestimates of output 
flux at low fluxes. The bias also increases at large off-axis angles. The best fit of the biases are 
shown as dotted lines for off-axis angle intervals 0'-2.5', 2.5'-4', 4'-6', 6'-8', and > 8'. 

Fig. 11. — Cumulative LogN-LogS for the soft and hard bands. The la error is shaded. Dash- 
dotted line represents the best fit from Moretti et al. (2003). Hard band LogN-LogS from Moretti 
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et al. is rescaled to that of 2 — 8 keV, assuming T = 1.4. 

Fig. 12. — Differential LogN-LogS for the soft and hard bands. The unit of dN/dS is number per 
10"^^ ergs cm~^ s"^ . Best-fit power laws are shown as solid lines. Dotted line represents a power 
law with a fixed index of —2.5. The data do not constrain the slope at high fluxes well. Dashed 
line shows the best fit of the hard band LogN-LogS from the SEXSI survey (Harrison et al. 2003) 
and the dash-dotted line is the best fit from Cowie et al. (2002). 

Fig. 13. — Hardness ratio vs. full band flux for the CLASXS sources. Open circles with arrows 

represent the upper or lower limits. Dashed lines with numbers label the hypothetical spectral 
indices, assuming the source spectra are single power laws with only Galactic absorption. Dotted 
line represents the typical error size of the hardness ratio for a source with hardness ratio of 1. 

Fig. 14. — Hardness ratio of the stacked sources in different flux bins. Crosses are the CLASXS 
sources and diamonds are the combined CDFs sources (Alexander et al. 2003). Sources with fluxes 
below 8 X 10~^^ ergscm"^ s~^ in the CLASXS catalog and 1 x 10~^^ ergscm"^ s~^ in the CDFs 
catalogs are not included to avoid incompleteness. Dashed lines are as in Figure 10. 

Fig. 15. — Light curves of the sources detected to be variable. The fluxes are normalized to the 
mean of all the observations. Numbers in the plots are the source numbers in the catalog, (a) Soft 
band; (b) hard band; (c) full band. 

Fig. 16. — (Upper panel) Fraction of sources that are variable in different flux bins. (Lower panel) 
Number of variable sources (dashed histogram) and total number of sources tested for variability 
(solid histogram) in the same flux bins as in the upper panel. 

Fig. 17. — Excess variability for the variable sources in each energy band vs. the flux of that band, 
(a) Soft band; (b) hard band; (c) full band. 

Fig. 18. — Spectral variability vs. full band fluxes for all the variable sources. The fluxes are in 
units of 10~^^ ergs cm~^ s~^ . Numbers on top of each plot are the source numbers in the catalog. 

Fig. 19. — Adaptively smoothed X-ray images of the extended sources superposed on R band 
images. (a) Sources 1 and 2; (b) Source 3; (c) Source 4. 

Fig. 20. — Regions for spectral extraction of Sources 1 and 2 on the Gaussian smoothed gray scale 
map of the clusters. The Gaussian kernal size is 6". Source regions are shown as circles. Elliptical 
annulus region is for the background extraction. 

Fig. 21. — X-ray spectra and best-fit MEKAL models of Source 1 (dash-dotted line) and Source 2 
(solid line). 

Fig. 22. — Combined probability contour of the temperature of Sources 1 and 2. Contour lines are 
1, 2, and 3a confidence levels. Cross is the best-fit temperature. Solid line represents the equality 
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of the temperature of the two clusters. 

Fig. 23. — X-ray spectrum of Source 3 and the best-fit MEKAL model. 

Fig. 24. — Curve-of-growths for the extended sources (Sources 1-4 shown in pannels (a) -(d)) 
normalized to the best-fit background. Dotted line shows the best fit of an integrated 2 dimensional 
Gaussian. 

Fig. 25. — Color-magnitude plot for the galaxies within 0.5' of the X-ray center. Solid lines show the 
model red sequence from Yee & Gladdcrs (2001) at the redshifts that best match the observations 
in Source 1 {z = 0.5), 2 {z = 0.5), and 4 (z = 0.45). In the plot for Source 4, the red sequences for 
z = 0.5 (lower solid line) and z = 1.0 (upper solid line) are shown. 

Fig. 26. — R band image of the gravitational lensing arc found associated with Source 3. 
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